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1 Introduction 

We have recently shown how spin-1/2 vectors and matrices can be derived from 
first principles [1, 2]. The method used is based on the Lande interpretation 
of quantum mechanics [3,6]. This essentially consists in using the probability 
amplitudes which characterize spin measurements between different directions 
as the basis set for the expansion of the same probability amplitudes. The 
f"*^ ' method is valid for all values of spin, and we have illustrated it by applying it to 

f^ . the case of spin 1/2. We have thereby obtained not only the Pauli spin matrices 

C) ' and vectors but the most general forms of these quantities as well. 

The subject of this paper is the derivation of the spin vectors and operators 
for the case of spin 1 . This task is far from routine because to the best of our 
knowledge, the generalized results we shall introduce are new. At the same 
time, we confirm the validity of the method. The fact that it works for both 
spin 1/2 and spin 1 indicates that it will apply to any value of J . Our success 
O^l ' in this program shows the credibility of the Lande interpretation of quantum 

^ . mechanics. 

In this paper, we demonstrate how to obtain the standard forms of the 
spin -1 operators and eigenvectors from first principles, and how to generalize 
\^5 ' their forms. After deriving the standard general formulas for the probability 

f^ . amplitudes that describe spin-1 measurements, we use them to find the most 

^N ' generalized forms of the operators for the a;, y and z components of spin and their 

eigenvectors. We also obtain the generalized ladder operators. The correctness 
f~i ' of these quantities is proved by the fact that in the appropriate limit, they yield 

Mh. the standard forms. The operators for the a;, y and z components satisfy the 

commutation relations, and the generalized ladder operators have the required 
H , effect on the eigenvectors of the generalized z component of spin. In addition 

I— { ' to these operators, we derive the operator for the square of the spin. As a final 

^^, step, we summarize the procedure to be used in applying our method to any 

value of J. 

Our method is justified fully and developed step by step in ref. 1. Here we 
shall use the method with a minimum of the underlying theory. 

This paper is structured as follows. In Section 2, we review those results of 
the Lande approach which we shall need to use. Here, we give the general form 
of the matrix eigenvalue equation in this approach, as well as the form of the 
matrix operator and its eigenvectors. In Section 3, we apply these results to the 
case of spin 1. Thus, in Subsection 3.1, we derive the forms of the operators and 
vectors for spin 1 in terms of implicit probability amplitudes for spin projection 
measurements. In Subsection 3.2, we summarize all the possible forms that the 
operators and vectors can take if we vary some of the parameters on which they 
depend. 

We show in Subsection 3.3 that we are able to derive the standard forms 
of the operators and vectors without explicit knowledge of the probability am- 
plitudes. We then derive the explicit forms of the probability amplitudes in 
Subsection 3.4, and verify their correctness in Subsection 3.5 by showing that 
they obey the Lande expansion for probability amplitudes, which the whole 
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theory is based on. We obtain the corresponding probabiUties from them in 
Subsection 3.6. Armed with the exphcit probabiUty ampUtudes, we obtain the 
generahzed form of the operator for z component of the spin as weU as its 
eigenvectors in Subsection 3.7. 

In Subsection 3.8, we give the generahzed forms of the x and y components 
of the spin and their eigenvectors. We observe that the three generahzed com- 
ponents obey the usual commutation relations. In Subsection 3.9 we give the 
ladder operators, which are found to perform their lowering and raising functions 
on the generalized eigenvectors of the z component of the spin operator. We 
derive the operator for the square of the spin in Subsection 3.10. In Subsection 
3.11, we list some properties of the generalized probability amplitudes. 

In Section 4, we summarize the steps to follow in order to apply this method 
to higher values of J. We conclude the paper after a discussion in Section 5. 

2 Basic Results 

Let three observables A, B and C belong to the same quantum system. Let each 
one have N eigenvalues. If the system is known to be in a state corresponding 
to the eigenvalue Ai of A when B is measured, then the probability amplitudes 
corresponding to the different values of B being obtained are xi^il^k) where 
A: = 1,2, ...,7V. But if C is measured, the probability amplitudes corresponding 
to the different values of C being obtained are ip{Ai; Ck) where k = 1,2, ...,N. 
Finally, if C is measured after the system has been ascertained to be in a state 
corresponding to Bi the probability amplitudes for obtaining the different values 
of C are 0(5,; Ck) where fc = 1, 2, ...,N. 

The probability amplitudes obey two-way symmetry: hence they satisfy the 
Hermiticity condition 

'iP{A;Ck)=nCk;A,). (1) 

They are related by the Lande expansion [3-6] 

N 

V-CA.sCfc) = ^x(^^;S,)'/'(B,;C:fc). (2) 

i=i 
Using this basic formula, we find that the matrix representation of the prob- 
ability amplitude ip{Ai; Ck) in the basis 4>{Bj; Ck) is [1] 



bP{A^■,Ck)] 



( x(A,;i?i) \ 
X{M;B2) 



(3) 



V X{A^■,BN) j 

This formula is essentially the one found in standard texts. It differs in giving 
the elements of the vector as probability amplitudes connecting an initial state 
corresponding to eigenvalue Ai with a set of intermediate states corresponding 
to another observable - in this case B. 



Consider now the differential eigenvalue equation 

A(x)V'(Afc; x) = AfcV(Afc; x) (4) 

where A^ is an eigenvalue, and x is a general quantity which represents the 
final eigenvalue, and whose spectrum is continuous. This eigenvalue equation 
emphasizes the idea that the eigenfunction is a probability amplitude connecting 
an initial state corresponding to the eigenvalue Afc with a set of final states whose 
eigenvalues are given by x. We use Eq. (0) to express ^{\k] x) as an expansion: 



N 



V'(Afc;x) = ^x(Afc;Sj)0(Bj;: 



(5) 



i=i 



Substituting this in the eigenvalue equation Eq. (|j), we obtain [1] the matrix 
eigenvalue equation which is the equivalent of the differential eigenvalue equation 



/ All ~ Afc Ai2 
A21 A22 — Afc 



V Ami 



A 



N2 



AiN \ 



i2N 



( x(Afc;Bi) \ 
x(Afc;B2) 



= 0, (6) 



ANN-\k j \ x{^k\BN) j 



where 



Anj = {HB,n;x) \A{x)\ (f>{Bj;x)) . 
The expectation value of the quantity R{x) is given by the formula 

N N 
z=l i=l 



where 



/ x(Afc;Si) \ 
x(Afc;S2) 



[V-fc] 



V x(Afc;Bjv) J 



and 



/ -Rii R12 
R21 R22 



[R] 



RlN \ 

R2N 



\ Rni Rn2 ■■■ Rnn j 
The matrix elements of [i?] arc 



(7) 



(8) 



(9) 



(10) 



(11) 



3 Specialization To Spin 1 Systems 

3.1 Derivation of General Matrix Mechanics Formulas 

The Landc formula Eq. (||) applies to any three sets of probability amplitudes 
which relate different observables of a quantum system. As we showed in ref. 
1, the formula is particularly useful for the treatment of spin. The reason is 
that spin projections with respect to different reference vectors have the status 
of different observables. Yet, the three sets of probability amplitudes have the 
same functional form. This allows us to deduce the probability amplitudes 
easily. 

Consider a system whose spin is known to be up with respect to the unit 
vector a. We are interested in the probability amplitudes for measurement of the 
spin along the new direction determined by the unit vector c. We are interested 
in the quantity R{cr ■ c), which is a function of the spin projection <t • c. We 
may also measure the spin projection with respect to the unit vector b, and the 
value of the R{(t ■ b). 

We shall measure the spin in units of h. Thus, the possible values of the spin 
projection are +1, and —1. The probability amplitudes will be labelled with 
the magnetic quantum numbers corresponding to the initial and final directions; 
the quantization direction for each quantum number will be given as a subscript 
of that quantum number. Thus for measurements from the initial direction a 
to the final direction c, the complete set of possible probability amplitudes is 

'0(m-j , rrif), with all possible nine combinations ni^ , rrif — —1, 0, 1. 

For measurement from the initial direction a to the final direction b, the 
probability amplitudes are x(fri^ ; rrif ), again with TOj , to , = —1, 0, 1. 

For measurement from the initial direction b to the final direction c, the 

probability amplitudes are (j>{m\ , rrif). In all cases, the ordering of the eigen- 
values is such that i,f = 1 corresponds to the eigenvalue +1, i,f = 2 corre- 
sponds to the eigenvalue and i, / = 3 to the eigenvalue —1. 

According to the Lande formula, Eq. (0), the following expansions hold for 
the probability amplitudes tp : 

+X(mf); (-1)(^))0((-1)(^);^^)), (12) 

Similar expressions hold for the probability amplitudes x s-nd 4'- The matrix 
representations of the probability amplitudes are 

mm^^;mf^)]^ x(mp;0(b)) (13) 



Let the quantity R{(t ■ c) have the value ri when the spin projection is 
+ 1, 7'2 when the spin projection is 0, and r^ when the spin projection is —1. 

Assuming that the initial state corresponds to the spin projection rnf'h,, the 
expectation value of R is 



(i?) = ^ ^(mrUi'^)) 



ri=l 



(14) 



In consequence of the expansions given by Eqn. (jlj) , the expectation value 
Eq. (|lj) may be written as 



where 



' X(mp^(+1)(^)) ^ 

Vx(™P;(-i)^^)) / 



mrnf'^)\ 



and 



R\\ Ryi Rl3 

[R] = ( R21 R22 R23 
R31 R32 R33 



The elements of [R] are 



(15) 



(16) 



(17) 



Ru 



<^((+l)(^);(+l)(^))| ri + |</,((+l)(b);0(^)) 
+ </.((+l)(^);(-l)(^))%3, 



r2 



(18) 



R12 = (/.*((+l)('');(+l)(^))^(O(b);(+l)(^))ri + +0*((+l)('');O(^))0(O('');O)(^))r2 
+0*((+l)(^); (-1)(^))0(O(^); (-l)(^))r3, (19) 

Ri3 = </.*((+l)(^M+l)'^^)0((-l)'^M+l)'^^)^i +0*((+l)'^^O(^))(/.((-l)(^);O(^))r2 
+</.*((+l)(^); (-l)(^)),/,((-l)(b); (-l)(^))r3, (20) 

i?2i = (/.*(O(^);(+l)(^')0((+l)(^);(+l)(^))ri + (/.*(o(^);O(^)),^((+l)(^';O(^')r2 
+0*(O(^); (-1)(^))0((+1)(^); (-l)(^))r3, (21) 

R22 = 0(O(S);(+1)(^)) %i+ 0(O(^);O(^)) %2+ (/.(0(^);(-l)(^') %3, (22) 



i?23 = (/.*(O(^);(+l)(^))0((-l)(^);(+l)(^))ri + (/.*(O('');O(^))</)((-l)('');O(^))r2 
+0*(O(^); (-1)(^))0((-1)(^); (-l)(^))r3, (23) 

R31 = 0*((-l)(^);(+l)(^))0((+l)(^';(+l)(^))ri +0*((-l)(^);oW)0((+l)(b);oW)r2 
+0*((-l)(^); (-l)(^))</)((+l)(^); (-l)(^))r3, (24) 



i?32 = r((-i)^''';(+i)'^^)0(o(^';(+i)^^^)^i + r((-i)^''';O(^')0(o('');O(^))r2 

+0*((-l)(^); (-l)(^))(/)(0(^); (-l)(^))r3 (25) 



and 



i?33= 0((-l)(''^(+l)('=)) ri+ </.((-!)(''); O^'^)) r2+ (/.((-!)(''); (-!)('=)) 
These matrix elements are given by 



(26) 
(27) 



3.2 Summary of the Various Possible Choices for the Ref- 
erence Vectors 

The forms that the quantities [tp] and [R] assume are determined by the choices 
we make for the vectors b and c. The various possibihties are Usted below. 

Case (a) : b 7^ a and c 7^ a. 

This is the most general case. In this case, the matrix representations are 



(a)._(c) 



[V'(mr;m;0] = 



X(mp;0(b)) 
Vx(mp;(-l)(b)) ] 



(28) 



while the elements of the operator [R] are given by Eqns. (|18| ) - (^ 
Case (b): b = a. 
The matrix representations are 



S,7li 1 



,(a)._(c)^ 



(29) 



while the elements of [R] are given by Eqns. (^^ - (|2^) with b replaced by a. 
Case (c): b = c. 



For this case, the matrix representations are 
while the operator is 



X(mp;0(^)) 



(30) 



[i?]= 

V 










Case (d): c = a. 

The matrix representations are 







(31) 



('^)._(=^)^ 



[ipiml ';my)] 



I x(mp;(+l)(^)) \ 

X(™p;0(^)) (32) 

Vx(mf^(-i)(^)) ; 

while the elements of [R] are given by Eqns. (|8|) - (|2^) with c replaced by a. 
Case (e): b = a and c = a. 
The matrix representations are 



[■0(TOj ; my )] 




(33) 



while the operator is 




(34) 



This case evidently corresponds to the standard forms of both the vectors 
and the operator, and the vector states happen to be the eigenvectors of the 
operator. 

3.3 Derivation of the Standard Forms of the Spin Vectors 
and Matrices 

We now employ the expressions for the elements of [R] to deduce the matrix 
forms of the spin operators. When R is the spin projection itself, then ri = 1, 
r2 = and r^, = —1. Cases (c) and (e) permit us to deduce the standard form 
of the operator for the spin projection along the axis of quantization. We shall 
call this the operator for the " z component" of spin, the quantization direction 
is not necessarily the z direction. We find that 



1 














-1 




[<yz\ =00 . (35) 

V -1 / 

We emphasize that the special form of " [az] " in Eq. ( p5| ) applies whenever 
we choose b = c, even if neither a nor b nor c coincides with the z direction. 
Nevertheless we find it natural to give this operator the subscript z. 

Now, since the eigenvalue of tr^ is 2, we have ri = r2 = rs = 2: the standard 
form of the operator for the square of the spin is given by Cases (c) and (e) as 



(36) 



For the two cases b = c and a = b = c, for which " [ayj\ " is given in Eq. (pq), 
we can deduce the matrix forms of the operators for the x and y components of 
spin by first obtaining the ladder operators [(t+\ and [u-] . We obtain the ladder 
operators by their actions on the vectors of [cr^] . To make a distinction between 
vector states and eigenvectors, we denote the latter by [^^], where i stands for 
+ , or — . We have 

[a+] [e+] = 0; [<J+] Ko] = A/2[e+]; [a+] [U = ^M (37) 

and 

[a_] [U = 0; [a_] [i,] = V2[S,_]- [a_] K+] = V2[U. (38) 

We thus derive the standard formulas 

V2 
[a+] = I V2 \ (39) 



and 








[a-] = I x/2 I . (40) 



From 



and 



we obtain 



V2 



[a.]^\{[a+] + [aA) (41) 



K] = -^([^+]-['^-]), (42) 



and 




(43) 



(44) 



Thus, we have derived from first principles the operators for spin 1, just as we 
derived the operators for spin 1/2 from first principles in ref. 1. Now, these spin 
matrices are clearly not the most general operators for the description of spin 
measurements. To obtain these general forms, we need the explicit expressions 
for the probability amplitudes so that we can substitute them in the formulas 
Eq. (§1) - (H). 

3.4 Explicit Expressions for the Probability Amplitudes 

The standard matrix operator for the spin is 







— i 





i 





— i 





i 






cr]=i +j^ 1 1 +k^ i -i \ ■ (45) 



We now take the dot product of this operator with the unit vector a. If the 
polar angles of a are {9, ip) , then in Cartesian coordinates, a = (sin 6 cos (p, sin 6 sin ip, cos ( 
A measurement of spin along the new direction defined by the unit vector a gives 
the values 0,±1. We have 



/ cos 6* ^sin6'e-*'^ 



72' 

-ijsin6'e*'^ -i^ sin6'e~''^ | . (46) 

\ ^sin6'e''^ -cos 6* 

The eigenvalues of this operator are +1, and —1 with the respective normal 
ized eigenvectors 



K+]= 72 ^in^ I (47) 




Ko]=| cose I (48) 

and 



-^ sin 9e"^ 

v2 




(49) 



The vectors Eqns. (^ - (^ are the eigenvectors resulting from solving the 
eigenvalue equation Eq. (|6|). Therefore, they are of the form x(a/c;-Bj), where 
Ofc is an eigenvalue of the matrix, and corresponds to the state that precedes 
measurement, while Bj is an eigenvalue corresponding to another state. In 
general B is not the quantity being measured. The quantity B is the spin 
projection along another direction which is defined by some unit vector d, which 
for the moment is unidentified. We can see that in the vectors [S,^], [£,q] and 
[C-]; the elements are probability amplitudes. Thus, for example, consider the 
elements of [^_|_]. These are probability amplitudes; they refer to measurements 
from the initial state defined by the spin projection +1 with respect to the 
vector a to the states characterized by the final projections +1, and —1 with 
respect to the unknown unit vector d. Denoting the probability amplitudes by 
the symbol X: we deduce that 

x((+l)(^);(+l)(^))=cos2^e-*^ (50) 

^ii+lf);0f))^l=sm9 (51) 



and 



^{{+l)i-);{-l)('^))^,in'-e^<P. (52) 



Similarly, we deduce from [£_q] and [^_] that 

x(O(«);(+l)(^)) = -i=sin0e-'^ (53) 

x(o(^);O)(2))=cos0, (54) 

x(0)(^);(-l)(^)) = -i=sin^e'^ (55) 

X((-1)(^);(+1)(^)) = -Sin2^e-^ (56) 

^((_l)(a).o(3)) = _Lsin0 (57) 

and 

x((-l)(^^(-l)(^)) = -cos2^e^'^ (58) 

These expressions are not necessarily the most general probability ampli- 
tudes because d may not be an arbitrary vector. In fact, since generalized 
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probability amplitudes must contain angles pertaining to two different direc- 
tions, we can see that the probability amplitudes Eqns. ( pO| ) - ( p8[ ) are not 
candidates for this role. However, in view of the Lande formula Eq. (0), we 
can use Eqns. (Q) - (JS^) to eliminate the unknown direction d and hence de- 
duce the general form of the probability amplitudes. Thus, we introduce a new 
unit vector c defined by the polar angles [9' , ip'). From the eigenvectors of the 
operator [<t • c] , we deduce that 

x((+l)(^);(+l)(^))-cos2|e-^', (59) 

^{{+lf);Oi^))^ lysine', (60) 

x{(+lf^;{-l)^'^^) = sm'^-e^^', (61) 

X(0(^); (+1)(^)) = --Lsin^'e-^'^', (62) 

^(0(^);0(d))^cos6l', (63) 

^(O(c);(_i)(3))^_Lsin0V^', (64) 

X((-1)(^^ (+1)(^)) = -sin2 ^-e~^^', (65) 

^((_l)(c).o(3))^_Lsin0' (66) 

and 

x{{-lf);(-lf)) = -co,'^-e^^'. (67) 

We can now use the Lande expansion Eq. (0) to find the general probability 

amplitudes ip{Tn\ ;mf ). Using the Hermiticity condition Eq. (w) to reverse 
arguments in the probability amplitudes Eqns. (bm - (pT)), and using the fact 
that all three sets of probability amplitudes ^, x and (j) have exactly the same 
form, we find that 



V'((+i)(^^ (+1)^^^) = x((+i)^"^; i+iY'^^Mi+i)'-''^] (+1)^^^) + x((+i)(^^ o('^))(/.(o('*); (+i)(^)) 

+x{{+if'>; (-i)(^))</.((-i)(^); (+i)(^)) 



-COS — e ^^ ^ ^ + sm - sm —e^^ ^ ' + -smOswYO 
Z ^ 

The other probability amplitudes are 



= cos2-cos2-e"*^'^"'^'^+sin2-sin2-e*('^-'^'' + -sin6lsin6l'. (68) 

2 2 2 2 2 ^ ' 
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■^((+1)(^);0('^)) = — f«in2^^;„fl'pi(ip-tp') „r.o«2_sin/3'p-»(v-v') 



sm — smW e ^^ ■»- / _ cos — sint/ e 



V2' 2 2 

+ sin 6* cos 6*'], (69) 



X((+1)('^^ (-!)(•=)) = cos2 ^ sin^ |e-'(^-^') + sin^ I cos2 |e'^(^-^') 

-i sin 6* sin 6*', (70) 

X(0(^);(+1)(^)) = -^[-sin6icos2-e-^('^-'^')+sin6lsin2-e'('^-'^') 
\/2 2 2 

+ cos 6* sin 61'], (71) 

^(0(a);0(^)) = isin6'sin6l'e"*('^~'^') + i sin6lsin6'V('^"'^') + cos6lcos6l', (72) 

X(0(^);(-1)(^)) = ^[-sin6isin2-e-*('^-'^')+sin6'cos2-e'('^-'^') 

-cos 6* sin 61'], (73) 



;^((_l)(a).(+i)(c)) ^ sin^ - cos2 -e-'(^-^') + cos2 - sin2 -e^(^-^') 

-i sin 61 sin 6*', (74) 

X((-1)^^^;0(^)) = ^[-sin2-sin6''e-*('^-'^')+cos2-sin6l'e'(^-^') 
V 2 2 2 

-sin 61 cos 61'] (75) 

and 



^((_l)(a) („i)(c)) ^ sin2 - sin2 -e-^(^-^') + cos2 - cos2 -e^(^-^') 
^w J y\ J J 2 2 2 2 

+ i sin 61 sin 6*'. (76) 

Once having obtained the generaUzed expressions Eqns. (pq) - d79), we are 
able to deduce the vector d which we used to obtain them. This vector is 
evidently k, the unit vector defining the z direction. Thus if we replace {9, ip) 
hy (9 = ip = 0), in the generalized expressions, we recover Eqns. ( pO|) - ( pq) . 
The only difference to note is that the probability amplitudes Eqns. (pq), (pTh 
and Eq. (pq) corresponding to the initial state —a are multiplied by the phase 
factor —1. 
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3.5 Verification of the General Expressions 

If the probability amplitudes Eqns. (Q) - (^^ are correct, they should internally 
satisfy the Lande formula Eq. (||). Let e be a new unit vector defined by the 
polar angles {9", if"). Then according to Eqns. (pi) - (vw , 



^((+l)(a); (+i)(e)) ^ cos^ I cos^ ^e-*(^-^") + sin2 | sin^ ^e'i^'^"^ 

+ - sin 6* sin 61", (77) 



X((+1)(^^ 0(^)) = ^[sin^ ^ sin0"e'(^-^") - cos^ ^ sin0"e"*(^-^") 
\/2 2 2 

+ sin 6* cos 6*"], (78) 



X((+1)(^);(-1)(^)) - cos2-sin2— e-'(^-^")+sin2-cos2— e*(^-^") 

-i sin 6* sin 61", (79) 



3" fl" 



Y(0(^^(+l)(^h = ^[-sin^cos^— e-'('^-'^")+sin6isin2— e'(^-v'") 
AV 7V y / ^2^ 2 2 



(80) 



^(0(a).o(«)) = isin6'sin6l"e-'('^-'^") + isin6'sin6l"e*('^-'^") 

+ cos 61 cos 61", (81) 



(0(^^(-l)(^)) = ^[-sin6lsin2^e-'('^~'^")+sin6'cos2^e'('^-'^") 
V 2 2 2 



-cos 61 sin 6*"], (82) 

X((-1)(^^; (+1)(^)) = sin2 6i/2cos2 Lg-'^'^-'^") + cos^ - sin^ Le*('^-<^") 

-i sin 6* sin 61", (83) 

X((-1)(^);0(^)) = ^hsin2-sin0"e^(^-^")+cos2-sin0"e-'(^-^") 

-sin 61 cos 6*" (84) 
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and 



+ - sin 6* sin 61". (85) 

Wc should be able to obtain all the probability amplitudes for measurements 
from a to c by using these expressions in the Landc expansion. Making use of 
the Hcrmiticity condition Eqn. ffl), wc find that 

x((+i)(^); (+i)(^)) = x((+i)^^; (+i)(^))x((+i)(^^ (+1)^)) + x((+i)^^^o(^))x(o(^); (+i)(^)) 

+X((+1)(^);(-1)(^))X((-1)(^^(+1)<^^) 

fl' ft" ft' ft" 

= cos^ - cos^ ^e-^(^'-^") + sin2 - sin^ ^e^(^'-^") 

2 2 2 2 

+ i sin 6*' sin 6*" (86) 

the same result which is found by directly by putting the appropriate angles in 
the relevant formula Eq. (pq). Thus, the Lande expansion is satisfied for this 

probability amplitude: in fact, it is satisfied for all others in the set xi'^t : inf) 
. Therefore, the probability amplitudes are correct and consistent. 

3.6 Expressions for the Probabilities 

The probabilities corresponding to the probability amplitudes Eqns. 
are 



P((+l)(^);(+l)(^)) = x((+i)^^';(+i)^^') 



= cos - cos h sm - sm — 

^22 22 

+ -sin6'sin6''cos(i^'- (^)]^ (87) 



P{{+\);0^^) = i[l - (cos6'cos6'' + sine'sin6''cos((p' - ip)f] 



P((+l)(«);(-l)(^)) = [cos^^sin^l+sin^^cos^l 

— sin 61 sin 6*' cos(93' -(/3)]2, (89) 

P[q(^)-q(^)) ^ [cos9cos9' + sm9smd'cos{ip' -<fi)]'^, (90) 

P(0'"'; (+1)<^)) = P((+1)("';0(^)), (91) 
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P(0("';(-1)('^') = P((+1)("';0('=)), (92) 

P{{-lf^; (+1)(^)) = P{{+lf^; {-lf'>), (93) 

P((-1)(");0(^') = P((+1)("^0(^)) (94) 

and ^ ^ ^ ^ 

P((-l)('^); (-1)(")) = P((+l)'"'; (+1)'"^). (95) 

The probabilities corresponding to a particular initial state should sum to unity. 
We verify that 

3 _ 

5^P(mf);mf)) = l, (z = 1,2,3). (96) 

n=l 

3.7 Explicit Formulas For The z Component Of The Spin 
Operator 



We are now in a position to use the defining expressions Eqns. ( |18D - (|26|), to 
deduce the expressions for the elements of the generalized " z component" of the 
spin operator. Since this operator is referred to the direction c, we now denote 
it by [cr--]. In the defining expressions, we have probability amplitudes of the 

form 4'{ml , rrv?'). Let the angles of b be {9, ip) and those of c be {9\ </?'). We 
find that 

(cr^)ii = Pii = cos cos 6*' + sin 6' sin 0'cos(<y9 — </?'), (97) 

(o'^)i2 = ^12 = —?= [^ sin 6 cos 6' + cos 9 sin 6' cos((y3 — ip') 
" V2 

+isin6''sin((^-(^')]j (98) 

('T^)l3 = (^.)22 = {(Tz)zi = 0, (99) 

(o'^)2i = ^^[— sin^cos^' + cos6'sin0'cos((/3 — (p') 

" \/2 

-isin6i'sin((^ -(/?')], (100) 

(a-)23 = (o--)i2, (101) 

(a^)32 - (a^)2i (102) 

and 

(a-)33 = -(^-)ll. (103) 
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When we set b = k in the expressions above we obtain the elements of the 
operator Eq. (^). The only difference to note is that due to the change of 
notation in Subsection 3.7, the angles that have to be set to zero are and ip. 

We now need the eigenvectors of the generalized operator. According to 
the matrix eigenvalue equation Eq. (0), the elements of the eigenvectors are 
probability amplitudes. Therefore, they have the forms of the generalized prob- 
ability amplitudes. When we examine Eqns. (pi) - (p5), we discover that the 
case for which a = c while b is arbitrary corresponds to the situation in which 
the operator has its most general form, while at the same time the vectors are 
eigenvectors of the operator. Hence the eigenvectors of the generalized operator 
correspond to the case where the initial state in each eigenvector belongs to the 
vector c, while the final states belong to the vector b. This is why in order to 
recover the special vectors Eqns. ( p7[ ) - (|49|), we set b = k. We keep in mind 
that b defines the intermediate observable. 

The elements of the generalized eigenvectors are probability amplitudes of 

Since they are referred to the arbitrary direction c, the 



(c). 



the form x("^ 

eigenvectors are now denoted by [C-^^'']i [S,^''] a-nd [4^ '']. They are 



'f) 



(+)i 



40)i 



[e 



(+)i- 



^((+l)(c).0(b)) 



cos2 ^ cos^ |e"*('^'-'^) + sin^ ^ sin^ ^e^'^v'-v) + i sin6l' sin 61 \ 
^[- cos^ 1^ sin 0e-*('^'-'^) + sin^ 1^ sin 0e'(^'-'^) + sin 6*' cos 61] 
cos2 ^ sin^ |e-*(^'-¥') + sin^ ^ cos^ le*^'^'-'^) - i sin 6*' sin 61 / 

(104) 



(0)l _ 



[e 



and 



/ x(O(^);(+l}(^0 
V x(o(^);(-i)(^)) 



cos2 I sin6''e~'('^'-'^) + sin^ ^ sm( 



y^iiv'-v) _|_ 5 



sm f cos ( 



i sin 9 sin O'e-^'^'^'-'^'^ + i sin 6* sin 9'e''-^'-f^ + cos 6 cos 6' 
^[-sin2|sin6''e-*('^'-'^) + cos^ | sinSV^'^'-'^) - sin g cos 61'] 

(105) 



sin2|cos2fe-^(^'-^) 



+ cos^ sm |e 



[^Ir^] = I x((-l)^^';0(^') 1 = 1 ^[-sin2|-sin6'e-*('^'-'^)+cos2Csin6'e^('^'~'^) 



sm^ ^ sin^ |e-»('^'-'^) + cos^ |- cos^ ^e'^v'-v) + i gin 6*' sin ( 



(106) 



x((-i)('=);(+i)(^)) 

^((_l)(c).0(b)) 

x((-i)(^); (-!)('')) 

Direct substitution verifies that these vectors do satisfy the eigenvalue equations 

K][4''^] = (+i)[4''^], (107) 

K] Kf ] = (108) 



^^sin^fe'^'^'-'^) - i sin 6*' sin ( 
2 e^ 

2 



smf cost 



and 
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3.8 Generalized Expressions for the x and y Components 
of the Spin 

We now seek the operators [<Jx] and [oy]- The most obvious way to obtain 
them is to use the ladder operators, which we can derive by their actions on 
the eigenvectors of [cr^]- This is the procedure foUowed in ref. 1 for spin 1/2. 
However, it proved time-consuming in that case, and would be even more so 
in this case. Another method would be to used rotations. However, neither of 
these two methods is necessary. Fortunately, as we have shown in ref. 2, we can 
get the operators [ax] and [uy] simply by changing arguments in the expression 
for the generalized operator [a^]. By the same change of arguments, the vectors 
of [cr^] transform to the vectors of [ux] or [uy]. We proved this procedure by 
applying it to the case of spin 1/2 for which the operators [ux] and [ay] were 
already known [2]. 

For the case of spin 1/2, we found that the transformation 9' ^ 9' — -k 12 
applied to [a^ and its eigenvectors yielded [a^ and its eigenvectors. The trans- 
formations 9' = 7r/2 and tp' —>■ tp' — tt/2 gave [ay] and its eigenvectors. When we 
try these transformations in this case, we find that we obtain the wrong results 
for both [ax] and [ay]. We know this because if we set 9 — 9' and ip = <p' in these 
operators in order to recover the standard results Eqns. (Esl) and (E5), we fail 
to do so. However a little trial and error shows that if we set 9' -^ 9 + 7r/2,we 
do obtain [ax] and its eigenvectors. Similarly, 9' = tt/2 and ip' ^' ip' + 7r/2 lead 
to the correct results for [ay] and its eigenvectors. The standard results are 
obtained in the limit 9 = 9 and p = p' . 

Using these transformations, the elements of [ax] are found to be 

{crx)ii = — cos sin 6*' -f sin 6' cos 6*' cos(i^ — p'), (HO) 



{ax)i2 = ^p[sin6'sin0' -I- cos6'cos6''cos(iy9 — (/?') 
v2 

+icoB9' sm{p — p)]^ (111) 



(fa;)2i = — ^[sin6'sin0' + cos6'cos6''cos((p — (/s') 

v2 

^% (112) 

(113) 
(114) 
(115) 





—i cos 9' sin((^ 


(cra;)23 = 


= (cra;)l2, 


(CTa;)32 = 


= (cr2;)21, 


(ct3;)33 = 


= -(CT:.)i1 



and 



{<yx)lZ = (CTx)22 = (0'x)31 = 0. (116) 

The eigenvector corresponding to the eigenvalue +1 is 
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[i. 



(+)l 



/ i(l - sine') cos2 |e-*('^'-¥') + i(l + sin0') sin^ §e^('^'-'^) + i si 



sm W cos ( 



^[-i(l - sin6'') sin6'e-'('^'-'^) + i(l + sin6'') sin6'e*('^'-'^) + cos6lcos6l'] 



V2 



\ i(l-sin6l')sin2|e-*('^'-'^) + i(l + sin6'') cos^ |e*('^'-'^) - isin6'cos( 

(117) 
The eigenvector corresponding to the eigenvalue is 




- cos 9' cos^ |e-^('^'-'^) + cos 6*' sin^ ^e'^f'-f) 



sm W sm 



0'] 



cos 9' sin 6'e~*('^'"'^) + i cos 6*' sin We 



i(v'-¥') 



The eigenvector corresponding to the eigenvalue —1 is 



^y cos6''sin2|e~*('^'-^) + cosei'cos^ f e*^*''-'^) + sine* sin 0'1 



(118) 



a 



(-)i 



i(l + sin6i')cos2|e-*(^'-'^) + i(l - sin6'') sin^ f e*(*''"^) - isin6lcos6l' 
1 [-1(1 + sine') sin 0e-*('^'-'^) + i(l - sin^') sin^f 



V2 



■ i(v'-v) _ 



cos f COS ( 



i(l + sine') sin^ |e-*(v''-^) + i(l - sine') cos^ |e*(*''-^) + \ sir 



(119) 
We confirm that these eigenvectors satisfy the eigenvalue equations 

M[d+^] = (+i)Ki+^], (120) 



['^M. 



(")i - 



and 



The elements of [ay] are found to be 

{o-y)ii = sin e sin((/7 - (/?'), 



(f^y)23 
(c^y)32 



1 

(0-^)22 = (o-?,)3i = 

1 

71 

(0'?/)l2 
(f^?/)21 



[cosesin((^ — (;9 ) — zcos(<y9 — (/?)], 

)22 = (o'?/)31 = 0, 

[cosesin((^ — If') + icos{ip — (/?')], 



and 



(CTy)33 = -(o-y)ll- 

The eigenvector corresponding to the eigenvalue +1 is 



(121) 
(122) 

(123) 

(124) 

(125) 

(126) 

(127) 
(128) 

(129) 
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/ ihicos2|e-*('^'-'^) +zsin2|e*(^'-'/') +sm6l] 
[4+^]= ^[§sin6'e-*('^'-'^) + |sin0e*('^'-'^)+cos6l] |. (130) 

\ i[-isin2fe-*('^'-'^)+icos2|e*(*''-'^) -sin 61] 

The eigenvector corresponding to the eigenvalue is 

^[zcos2|e-*(^'-^)+zsin2fe*(^'-^)] 
K^] = I ^[-ismde-'(v'-'P) +ism9e'(v'-v)] |. (131) 

^[isin^ fe-*('^'-'^) + icos2 |e*('^'-'^)] 

The eigenvector corresponding to the eigenvalue —1 is 

i hi cos2 |e-*('^'-¥') + i sin^ |e'('/''-¥') - sin ( 
[4"'] = I ^[isin6'e-^('^'"'^) + |sin6'e*('^'-'^)-cos0] |. (132) 

i[-isin^ «e-*('^'-'^) + icos^ fe^(^'-¥') + sin( 

Direct substitution verifies that these eigenvectors satisfy the equations 

Wymi;^] = i+mi'\ (133) 



and 



K]Kf]=0 (134) 

Wymi-^] = {-mi-^ (135) 

Calculation shows that the generalized operators [cfx], [<7y] and [ctz] satisfy 
the commutation relations 



[[(^x], [t^]] = i[(Jz]; [[(Jy], [cr-]] = i[(7x] and [[a^], [a^]] = i[ay] (136) 

and that the sum of their squares gives the matrix for the square of the total 
spin: 



krl^ + \cr„f + W 




(137) 



3.9 The Ladder Operators 

The elements of the raising operator [(j+\ are 

((j+)ii = — cos6'sin6'' + sin6'cos6''cos(</? — (ys') + zsin6'sin(</? — </?') (138) 

(o'+)i2 = —p [sin 6* sin 0' + cos 6* cos 6*' cos((p — (^') + cos((^ — </?') + zcos6''sin((^ — ip') 
v2 

+icos6'sin((/3-(^')] (139) 
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(cr+)2i = —=[sm6sin9' + cos cos 0' cos((p — (p') — cos{ip — ip') — icos9' sin{ip — ip') 
+icos9sm{Lp-'p')] (140) 

(cr+)23 = (ct+)i2 (141) 

(^+)32 = (ct + )21 (142) 

(cr+)33 — COS 6* sin 0' — sm 6 cos 6' cos{(p ~ ip') — isin6'sin(<^ — ip') (143) 

and 

(ff+)l3 = (^+)22 = (ct+)31 = 0. (144) 

The elements of the lowering operator [ct-] are 



(cr_)ii = — cos6'sin6'' + sin6'cos6''cos((/3 — (/?') — isin6'sin((^ — (^') (145) 

(ct_)i2 = — ^ [sin sin 0' + cos cos 0' cos ((/j — ip') — cos{ip ~ ip') + icosO' sm{ip — ip') 
v2 

—icos9sin{p — p')] (146) 



(cr_ )2i = ^p [sin 6* sin 6* +cos6'cos6' cos(99 — </?')+ cos((/7 — </?') — i cos sin(<y9 — </?') 
v2 

-icos6'sin((^ -(/?')] (147) 

(fT-)23 = (a-)i2 (148) 

(fT-)32 = (fT-)21 (149) 

(cr_)33 = COS sin 6*' — sin 6* COS 6*' COS ((y9 — p') + zsin0sin((/? — p') (150) 

and 

(f7_)i3 = (a_)22 = (a-)31 = 0. (151) 

We confirm that when these ladder operators act on the generahzed vectors 
of [a^],Eqns.^-(|l|), they give 

[a+][4+)]=0; [a+][ef] = V2[4+^]; [a+mt^] = V2[^i°^] (152) 

and 

[a_][ei-'] = 0; [a_][ei°)] = V2[ei-^]; [a_][ei+)] = V2[ci"^ (153) 
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3.10 Generalized Form for the Square of the Spin 

We have seen that we do indeed obtain the elements of the matrix for the square 
of the spin by the laborious procedure of squaring its components and adding. 
Actually, this labour is unnecessary. We can directly infer that the generalized 
form of this operator is diagonal. 

The elements of the matrix for the square of the spin are given by the general 
expressions Eqns. ( p^ ) - (26). Now, the square of the spin has the same value 
no matter what the spin projection is. For such a quantity, the matrix operator 
is diagonal. We shall now prove that the off-diagonal elements all vanish by 
showing this for Ri2. 

According to Eqn. (E9), 



+0*((+l)(^); (-1)(^))0(O(^); i-lf^)r3 (154) 

Using the Hcrmiticity condition Eqn. ffl), and the fact that ri = r2 = r^ = r, 
we see that 

i?i2 = [(/)(O(^);(+1)(^))0((+1)(^);(+1)(^)) + 0(O(^);O(^))(/)(O(^);(+1)(^))) 
+0(O(^); (~1)(^))0((-1)(^); {+lY%. (155) 

Thus, according to the Lande expansion Eqn. (g), Eqn. (155) is 

i?i2 = 0(O(^);(+l)(^))r = O. (156) 

The same reasoning shows that all the other off-diagonal elements vanish. 
On the other hand, if we neglect the eigenvalue which multiplies each one, the 
diagonal elements are just the sums of all the probabilities corresponding to the 
same initial state. Thus for example. 



Rii 



0((+l)('');(+l)(^)) + 0((+l)('');O(^)) + 0((+l)('');(-l)(^)) 



(157) 

Each diagonal element therefore equals the eigenvalue. In consequence, we find 
that 




(158) 



a result we obtained earlier by squaring the components of the spin and adding. 
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3.11 Further properties of the probability amphtudes 



Apart from the Hermiticity condition, the probabiHty amplitudes satisfy the 
following relations. Given two reference directions defined by the unit vectors 
a and c, we readily verify that 



X*((+l)(^); (+1)(^)) - x((-l)(^); (-1)^^'), 



(159) 



and 



^*((+l)(a).o(c) 

X*((+l)(^);(-l)(^' 
^*(0(a);(+i)(c) 

^*(0(«);0(^) 
^*(0(a).(_l)(c) 

X*((-l)(^);(+l)(^' 

^*((_l)(a).0)(c) 



= -x((-l);,0('^)), 

= x((-i)(^);,(+i)(^)), 

= -x(o(^^,(-i)(^'), 

= x(o(^^o(^)), 

= -x(0(^';(+l)(^)), 

= -x((+i)^^^;o(^)) 



x*((-i)(^); (-i)(^') = x((+i)^^'; (+1)^^^). 



(160) 
(161) 
(162) 
(163) 
(164) 
(165) 
(166) 

(167) 



4 Results for Arbitrary Values of J 

The method we have presented of obtaining generalized spin vectors and oper- 
ators has now been applied to the cases of spin 1/2 and spin 1. It is useful to 
summarize the steps to follow in obtaining the generalized results for any value 
of J. 

Given any value of J, the first order of business is to label the 2 J+ 1 projec- 
tions so that, starting from the maximum projection, the index 1 corresponds 
to the projection Jh, the index 2 to the projection (J — l)h, and so on until 
the index 2 J + 1 which corresponds to the projection —Jh. Next we obtain the 
standard form of the operator [cfz]- In units of h, this will be 



f J 
J-1 





V 






-J + 1 





\ 




-J J 



(168) 



The normalized vectors of this operator will be 



[C. 



( M 



[o ) 



K,7-l] 



1 
[o ) 



K-,/] 



/ \ 



['l ) 



etc. 



(169) 
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corresponding to the eigenvalues J, (J — 1),....,~J respectively. 

Next we deduce the ladder operators by means of their action on the eigen- 
vectors of the z component. From them we get the x and y components of the 
spin. Using the three components of the spin, we obtain the operator 

k • a] = [a--] = (i[o-^.] +j[(Ty] + k[cr;j]) • a 

where a — (sin 9 cos ip, sin sin ip, cos 9) . 

The eigenvectors of this operator will contain as their elements probabil- 
ity amplitudes of the form xi^Tii l^f ), where d is a vector which has to be 
deduced. Similarly, if we now seek the eigenvectors of the operator [o" • c] , where 
c = {sm9' cos (p', sin 9' sin ip' , cos 9'), we obtain all the probability amplitudes 

x(m^'^ ; nif ) because they will be the elements of the eigenvectors of this op- 
erator. We then use the Lande expansion Eqn. (0), to eliminate d in order to 
obtain the generalized probability amplitudes which represent spin projection 
measurements from a to c. 

Armed with the generalized probability amplitudes, we obtain the general- 
ized operator for the z component of the spin. We have to keep in mind that 
these are defined in terms of an intermediate reference vector b, and the final 
reference direction c with respect to which we seek the spin projection. Let the 
spin projections with respect to the vector b be labelled from the largest to the 
smallest as 6i, &2v-^2j-i-i- Let the spin projections with respect to the vector c 
be labelled in the same way. Then the elements of [a^] are 

2,7+1 

('^c)^ = X! (t){hf,Ck)(t){ck;bi)rk. 
fc=i 

In this formula, r^ is the fcth projection of the spin with respect to the 
direction c, measured in units of ft. 

The eigenvectors of the generalized operator [a^] contain as their elements 
generalized probability amplitudes. In order for them to be eigenvectors of [cr--], 
the initial direction in the probability amplitudes must be c, while the final must 
be b. We arrange the amplitudes into column vectors with all those starting 
from one projection belonging together. With regard to the final direction, 
the probability amplitudes are arranged so that the one corresponding to the 
maximum value of the projection comes first, then the one corresponding to the 
next largest projection and so on. 

With both [cr-] and its eigenvectors known, we obtain [ux] and its eigenvec- 
tors by applying the transformation of the kind Lp' ^> Lp' + 1^/2 or ip' ^ Lp' — 7r/2 
on these quantities. It would appear from the two cases so far seen that in 
general, a little trial and error might be necessary. We have selected the correct 
transformation if when we set 9 = 9 and ip = ip' we obtain the standard results. 

Transformations of the form 9 = 7r/2 and ip = ip' + 7r/2 or ip = ip' — 7r/2 
will give us [<7y] and its eigenvectors; again the requirement that the standard 
forms be obtained in the limit 9 = 9 and ip = ip' will determine which of these 
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transformations is correct. With all three components known, the raising and 
lowering operators can also be obtained. 

5 Discussion 

In this paper we have extended the ideas we first introduced in ref. 1 by using 
them to derive generalized operators and vectors as well as probability ampli- 
tudes for spin-1 systems. That we have succeeded in doing so lends strong 
support to the ideas of Lande, which underlie our work. We consider it quite re- 
markable and important that the matrix treatment of spin can be derived from 
a probability amplitude basis. It should be remembered that spin is normally 
represented as being fundamentally of matrix nature. But here we have given 
spin a treatment analogous to that given to orbital angular momentum whose 
matrix description is based on the spherical harmonics. This shows that in the 
treatment of spin, matrices can be dispensed with if necessary. Perhaps other 
quantities in quantum mechanics that are treated by a purely matrix approach 
can also be reduced to a probability amplitude basis. The case of isotopic spin 
comes to mind. 

Our method depends crucially on the idea that a probability amplitude must 
always contain two state labels, one label being for the state that obtains before 
measurement and the other for the state that results from measurement. If 
we accept that the main difference between the spin probability amplitude and 
the ordinary wave function is that the former refers to discrete final states 
while the latter refers to continuous final states, then we may interpret the 
wave function differently from the standard way. Since the wave function in 
coordinate space results as the solution of an energy eigenvalue equation which 
is cast in terms of the coordinates, it is seen to be the probability amplitude 
for obtaining various values of position if the measurement starts from a state 
characterized by the corresponding energy eigenvalue. Likewise, for a given 
energy eigenvalue, the wave function in momentum space is the probability 
amplitude for obtaining various values of the momentum upon measurement if 
the initial state is characterized by the energy eigenvalue. For a system with a 
spherically symmetric potential, the radial solution is the probability amplitude 
for obtaining various values of the radial coordinate if the measurement starts 
from a state defined by the energy eigenvalue. The spherical harmonics are 
probability amplitudes for obtaining various angular positions starting from a 
state characterized by the angular momentum corresponding to the values of I 
and m. 

We notice that the results we have presented in this paper also apply to I = 1 
orbital angular momentum. They obtain if we measure the angular momentum 
projection twice. This kind of measurement is seen to be different from that 
which the spherical harmonic represents. As stated already, a spherical har- 
monic connects an angular momentum state and a position state, while the 
probability amplitudes presented in this paper connect two angular momentum 
states or two spin states. 
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It is the practice to give the standard forms of the spin vectors when con- 
structing a complete wave function for a quantum system. We now reahze that 
in order for the Umited standard forms to be useful, we must be careful to use the 
generalized operators when actually computing such quantities as expectation 
values. We must match spin vectors and operators properly in order to obtain 
the correct results. The intermediate reference direction connects the quantities 
in the expressions for the vectors with the quantities in the expressions for the 
operator. The final direction in the expressions for the vectors must correspond 
to the initial direction in the expressions for the elements of the operator. In 
this connection, we should distinguish clearly between the vectors for general 
calculations of the expectation value and those which are eigenvectors of the 
operator. The former refer to an initial reference direction which is different 
from the final direction with respect to which we seek the spin projection. The 
latter on the contrary are such that the initial direction equals the final direc- 
tion. Thus the vectors that we must use if we seek to calculate the expectation 
value when the spin measurement is a repeat measurement are the eigenvectors 
of the operator. 

We have seen that matrix mechanics vectors do not directly give the proba- 
bilities of finding eigenvalues of interest in individual measurements. According 
to the Lande view, this is however what the wave function gives. Neverthe- 
less, the vectors are extremely compact because each vector does contain all the 
probability amplitudes. 

In conclusion, we have obtained from first principles new generalized formu- 
las for spin 1 operators and vector states by means of the Lande interpretation 
of quantum mechanics. It is our belief that a re-examination of other standard 
quantum results in the light of this approach will yield useful generalizations of 
these results and provide valuable insights into quantum theory. 
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